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ABSTRACT 

Some Lyman continuum photons are likely to escape from most galaxies, and these 
can play an important role in ionizing gas around and between galaxies, including gas 
that gives rise to Lyman alpha absorption. Thus the gas surrounding galaxies and 
in the intergalactic medium will be exposed to varying amounts of ionizing radiation 
depending upon the distances, orientations, and luminosities of any nearby galaxies. 
The ionizing background can be recalculated at any point within a simulation by 
adding the flux from the galaxies to a uniform quasar contribution. Normal galaxies are 
found to almost always make some contribution to the ionizing background radiation 
at redshift zero, as seen by absorbers and at random points in space. Assuming that 
~ 2 percent of ionizing photons escape from a galaxy like the Milky Way, we find that 
normal galaxies make a contribution of at least 30 to 40 percent of the assumed quasar 
background. Lyman alpha absorbers with a wide range of neutral column densities are 
found to be exposed to a wide range of ionization rates, although the distribution of 
photoionization rates for absorbers is found to be strongly peaked. On average, less 
highly ionized absorbers are found to arise farther from luminous galaxies, while local 
fluctuations in the ionization rate are seen around galaxies having a wide range of 
properties. 

Key words: diffuse radiation - quasars: absorption lines - intergalactic medium - 
galaxies: structure. 



1 INTRODUCTION 

The extragalactic background of Lyman continuum photons 
plays an important role in ionizing the many absorption 
line systems seen shortward of Lya emission in quasar spec- 
tra, including those seen at low redshifts using the ultra- 
violet capabilities of the Hubble Space Telescope (Bahcall 
et al. 1996). The ionizing background intensity determines 
the neutral gas fraction in Lya forest absorbers and the ion 
ratios seen in metal absorption line systems. Furthermore, 
understanding the ionizing background is important for de- 
veloping theories for the formation and evolution of galax- 
ies and the Lya forest and possibly for finding the bary- 
onic mass contained within galaxies and the intergalactic 
medium. Interesting questions remain as to how much of 
this background is contributed by galaxies and what role 
galaxies play in ionizing gas around and between them that 
is detected as Lya absorption. 

The only largely neutral Lya absorbers are the damped 
systems which have neutral hydrogen column densities 



Nhi > 10^" '' cm~^. Ionized absorbers include some Ly- 
man limit systems {Nhi > 10^^ '^ cm~^) and Lya forest ab- 
sorbers which have lower neutral column densities. Lyman 
limit systems are thought to arise around galaxies (Bergeron 
& Boisse 1991; Steidel 1995), while Lya forest absorbers 
have been detected which are as weak as Nhi ~ lO'^^ cm^^. 
Although many of these may be associated with the smallest 
amounts of intergalactic gas, including that in void regions 
(Dave et al. 1999; Stocke et al. 1995; Shull et al. 1996; Pen- 
ton, Stocke, & Shull 2002), at least some stronger forest ab- 
sorbers are found near luminous galaxies. In particular Lya 
absorption is almost always detected at a similar redshift to 
a galaxy which is found within ~ 200 kpc of a quasar line of 
sight (Bowen, Blades, & Pettini 1996; Lanzetta et al. 1995a; 
Le Brun, Bergeron, & Boisse 1996; Chen et al. 1998; 2001). 
Bowen, Pettini, & Blades (2002) have shown recently that 
while nearby Lya absorbers are difficult to match with par- 
ticular observed galaxies, Nhi is correlated with the local 
density of detected luminous galaxies. 

In carefully observed spiral galaxy discs the neutral col- 
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umn density falls off slowly with radius over most of the ex- 
tents as seen with 21 cm HI observations, but then a rapid 
truncation is seen. Bochkarev & Sunyaev (1977) first sug- 
gested that a truncation would occur in a spiral disc at a 
sufficient radius where the HI column becomes ionized by 
an extragalactic background. Disk edges have been modelled 
also by Maloney (1993), Dove & ShuU (1994a), and Cor- 
belli & Salpeter (1993). More recently ionized gas has been 
detected in Ha emission using a Fabry-Perot 'staring tech- 
nique' (Bland-Hawthorn et al. 1994) beyond the HI edges of 
several nearby galaxies (Bland-Hawthorn, Freeman & Quinn 
1997; Bland-Hawthorn 1998). 

The intensity of the ionizing background radiation has 
been measured using the 'proximity effect', or the paucity of 
absorption lines close to a quasar emission redshift, at low 
redshift first by Kulkarni & Fall (1993) and more recently 
by Scott et al. (2002) who find J(912A) = 7.613.0 x 10"^^ 
erg cm~^ s~^ Hz~^ sr~^, or a frequency and direction av- 
eraged ionization rate of 1.9 x 10~^^ s~^ at redshift z < 1. 
Some evidence for redshift evolution in the background is 
seen. Other methods for estimating this intensity give gener- 
ally consistent results typically within the lower end of their 
uncertainty range, including limits on Ha emission from 
high-latitude galactic clouds (Vogel 1995; Tufte, Reynolds, & 
Haffner 1998; Vogel et al. 2002) and extragalactic HI clouds 
(Stocke et al. 1991; Donahue, Aldering & Stocke 1995) and 
estimates from the HI galaxy disc edges (Maloney 1993; 
Dove & ShuU 1994a; CorbelU & Salpeter 1993). 

Using quasar spectra and considering reprocessing of 
photons by the intergalactic medium, Haardt & Madau 
(1996), Dave et al. (1999), and ShuU et al. (1999) have calcu- 
lated the history of the intensity of the ionizing background 
down to low redshifts. Their values are approximately con- 
sistent with the measurements above, so that quasars make 
an important and possibly dominant contribution to the ion- 
izing background even at low redshifts. While the number 
density of quasars is low at low redshifts, the universe be- 
comes optically thin to ultraviolet photons at redshifts z <^2 
(Haardt & Madau 1996) such that ultraviolet photons emit- 
ted at 2 ~ 2 are likely to survive without reprocessing until 
2^0. In contrast, the ionizing spectrum at higher redshifts 
is modified by absorption and reemission (Fardal, Giroux, 
& ShuU 1998). 

What contribution might galaxies make to the ioniz- 
ing background at low redshifts? Some suggestions have 
been made that insufficient numbers of ionizing photons 
(less than one percent) escape from galaxies for an impor- 
tant contribution to the ionizing background (Deharveng 
et al. 1997; Henry 2002), though others measure (Bland- 
Hawthorn & Maloney 1999; Leitherer et al. 1995; Goldader 
et al. 2002; Hurwitz, Jelinsky, & Van Dyke Dixon 1997) 
and model (Dove & ShuU 1994b; Dove, ShuU & Ferrara 
2000) higher escape fractions between three and ten percent. 
Giallongo, Fontana, & Madau (1997), ShuU et al. (1999) 
and Bianchi, Cristiani, & Kim (2001) find that star-forming 
galaxies could even dominate the ionizing background if at 
least a few percent of the ultraviolet photons escape. Some 
ionizing photons are likely to escape from most galaxies. 
Bland-Hawthorn (1998) suggests that the gas detected be- 
yond the HI disc edge in several spiral galaxies is ionized 
by stellar populations within the galaxies, as the emission 



measures are stronger than those predicted (Maloney 1993; 
Dove & ShuU 1994a) for an extragalactic background. 

Given that many stronger Lya absorbers are found close 
to galaxies, it is possible that stellar populations within 
galaxies make some contribution to the ultraviolet photons 
that ionize any nearby absorbers. Thus what is measured as 
ionizing background radiation may vary in intensity with lo- 
cation, depending upon the galaxy clustering environment or 
the properties, such as brightness and extinction behaviour, 
of any nearby galaxies. The method for simulating a fluc- 
tuating ionizing background is described in Section 2, while 
the resulting fluctuations are discussed in Section 3. Results 
of varying model parameters are discussed in Section 4, and 
Section 5 describes the relationship of the ionization rate 
fluctuations to the properties and locations of galaxies. The 
value of Hq is assumed to be 100 km s~^ Mpc"'^. 



2 METHOD AND SIMULATIONS 

The simulation used here is an updated version of that first 
described in Linder (1998), where in each case here 12590 
clustered galaxies are placed in a cube with an edge of 28.9 
Mpc (except where this edge is adjusted at redshift one). 
Lya absorbers arise in gas within and extending from galaxy 
discs, which are modelled as in Charlton, Salpeter, & Hogan 
(1993) and Charlton, Salpeter & Linder (1994). Each galaxy 
has an exponential inner disc and an outer extension where 
the column density declines as a power law with radius, as- 
sumed as galaxy discs are generally found to be exponential 
while absorbers roughly obey a power-law column density 
distribution at lower column densities (see Appendix^. In 
reality a smoother transition probably occurs, and some ev- 
idence for such a transition has been seen by Hoffman et al. 
(1993). The radius at which each HI disc changes from ex- 
ponential to power law decline was defined previously (Lin- 
der 1998) in terms of the disc ionization edge, as little is 
known about this switching radius. Since the ionizing back- 
ground varies here, however, it makes more sense to define 
this switching radius in terms of the galaxy disc scale length. 
Assuming the switch from exponential to power law occurs 
in each galaxy at a radius of four HI disc scale lengths, sim- 
ilar results are seen in the absorber counts arising as com- 
pared to absorber counts simulated using the previously de- 
fined switching radius. 

The galaxies are chosen to have visible properties based 
upon observed distributions of galaxy parameters, such as 
a Schechter luminosity function and a flat surface bright- 
ness distribution (McGaugh 1996). The HI disc scale lengths 
{fiHi) are assumed to be proportional to the B scale lengths 
(/i_g), where hni ~ 1.7/is to start, as this was previously 
found to give rise to reasonable absorber counts in Lin- 
der (1998). The gaseous properties of the galaxies can thus 
be related to the visible properties. This makes sense as 
galaxy discs are generally somewhat larger in HI than in 
optical images, although the HI sizes of galaxy discs also 
vary depending upon the location within a cluster (Cayatte 
et al. 1994) so that the average ratio of hs/hHi is uncer- 
tain. Each galaxy is placed in a cube of space, where the 
positions are chosen to be clustered as described in Linder 
(2000) using a fractal type method based upon Soneira & 
Peebles (1978). Random lines of sight through the box can 



Galaxies as Fluctuations in the Ionizing Background Radiation 3 



then be simulated, and an absorption line is assumed to arise 
when these lines of sight intersect any disc or outer exten- 
sion. Neutral column densities for absorbers are found by 
integrating the HI density along the line of sight. 

Each galaxy disc has an ionization edge, or radius be- 
yond which no layer of neutral gas remains. The vertical ion- 
ization structure of the gas is modelled as in Linder (1998) 
which is similar to the model in Maloney (1993). Inside of 
this ionization radius, the gas is assumed to have a sandwich 
structure, where the inner shielded layer remains neutral 
and has a height (zi) determined by equation (6) in Linder 
(1998). The gas above height (zi) and beyond the ioniza- 
tion radius is assumed to be in ionization equilibrium where 
Orecntot = C'^Hi & highly ionized hydrogen gas with total 
(neutral plus ionized) density ntot and neutral density uh- 
The recombination coefficient is a^ec = 2.42 x 10~^^ cm"^ 
s"^(r/10, OOOK)""'^^ and the gas temperature T is assumed 
to be 20,000 K. The frequency- and direction-averaged ion- 
ization rate is determined at each point in space while 
converging numerically upon the ionization edge (the mini- 
mum disc radius where {zi = 0)) in a self-consistent manner, 
and is calculated as described below. 

The intensity of ionizing radiation is allowed to vary 
at each point in space within the box. The ionization rate 
^, which was previously assumed to be constant in Linder 
(1998; 2000), is recalculated here at each point in space, 
where — + C,gai and C^i, is the contribution from quasars. 
The value of C,gai can be recalculated at any point within 
the box based upon the flux from the surrounding galaxies. 
The value of C,b = 3.035 x 10"^'* s~^ is assumed from the 
calculation of Dave et al. (1999) at 2: = based upon spectra 
from Haardt & Madau (1996). 

Bland-Hawthorn & Maloney (1999) modelled the escape 
of ionizing photons from our galaxy by extrapolating from 
a calculation of the ionizing photon surface density at the 
Solar Circle using nearby O stars by Vacca, Garmany, & 
ShuU (1996). Bland-Hawthorn (1998) gives a simple estimate 
of the number of ionizing photons escaping from the Galaxy, 



where Lp — 2 x 10 e '^r 
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phot cm s 



in the first equation of their Appendix, at some point with 
distance r, where Q is the angle from the galactic pole, and 
r is the Lyman limit optical depth. In order to add such a 
radiation field to ( in the units above, we need a frequency- 
averaged quantity which also considers the ionization cross 
section for hydrogen, where an ~ 6.3 x 10~^^ (v/uo)"'^ at 
frequency f, and vq is the Lyman limit frequency. Thus, 
for example, the contribution from our galaxy to at some 
point at a distance r^pc from its center would be 



Cga 



(1) 



where tp — if^dv. The mean value for an is weighted by 
assuming that ipi, oc v'"" , so that 

C,gal,MW = <^0-g(i/o) { ^^ j^^q) ' (2) 

We assume as = 2.5, which gives {(Jh) ~ (JH(vo)/'i, while 
Sutherland & ShuU (1999) prefer a similar Qs ~ 1.9 to 2.2 for 
starburst galaxies, which would result in {(Jh) ~ gh{i'o)/'^- 
Bland-Hawthorn & Maloney (1999) assumed that our 
galaxy has an axisymmetric exponential disc where the ion- 
izing photon surface density nd{r) = noe"'"'''' for disc scale 



length h and central ionizing photon surface density no. In- 
tegrating nd{r) over the disc area out to an infinite radius 
gives a number of ionizing photons which is proportional 
to , while the ionizing photon surface density can also 
be expressed as an ionizing surface brightness (/x) where 
n/no — lOO'-'^""''-'^^. Thus we extrapolate the formula above 
to other galaxies by correcting for variations in central sur- 
face brightness (/^o) and disc scale length for galaxy i as 
compared to Galactic (MW) values, by assuming that the 
ionizing scale lengths and central surface brightnesses are 
proportional to the B values. The galaxy contribution to 
for some galaxy i becomes 



C,gal.i — '-P 



(x) 



10 



(3) 



At any point in space at which we wish to calculate C, we sum 
the contributions from all the galaxies in the simulation, so 
that 



EC 



(4) 



where galaxy i is at a distance r, and Qi is calculated for 
each galaxy as in Appendix IHl 

When calculating the ionizing intensity seen by the 
outer part of a galaxy, extinction from gas and dust is im- 
portant in shielding absorbing gas from being ionized by the 
inner parts of the galaxy. Although some neutral gas must 
remain around galaxies which is often seen to give rise to 
absorption, the galaxy itself still may be the most impor- 
tant contributor to ionizing the gas in its outer parts. To 
start we assume that each galaxy disc is flat within two HI 
scale lengths and then warped by ten degrees beyond that 
value (Briggs 1990), so that the outer parts of the disc are 
exposed to some ionizing radiation which escapes from the 
inner regions. 

A value of r = 2.8 is preferred by Bland-Hawthorn 
(1998) when modelling our Galaxy, although the preferred 
value could be different as a result of an error (Bland- 
Hawthorn & Maloney 2001). On the other hand it has been 
difficult to detect any dust extinction in low surface bright- 
ness (LSB) galaxies (for example O'Neil, Bothun & Impey 
1997). Little is known about the gas-to-dust ratio in LSB 
galaxies, although the total extinction is likely to be higher 
in high surface brightness (HSB) galaxies. Thus it might 
make sense to assume that r,; is related to central surface 
brightness for the simulated galaxies. Thus a linear rela- 
tionship was assumed where r = 2.8 for a galaxy with a 
Freeman surface brightness value of /^o,i ~ 21.65, and r = 
for Ho A = 25, so that n = — 0.836(/^o,i — 25). 

Shadowing, or extinction from gas between an absorber 
and a given ionizing source, is not taken into account as this 
would require substantially more computing time. This is 
a reasonable approximation in the sense that the universe 
is optically thin at redshifts ^ 2 (Haardt & Madau 1996) 
Shadowing could make the calculated ( values slightly lower 
in some cases, although typically the second closest galaxy to 
an absorber contributes only a few percent of the ultraviolet 
flux that the closest galaxy contributes. 

Starburst galaxies are not treated as emitting differ- 
ently from other galaxies within the simulation, although 
the galaxy population is chosen within the simulation to be 
consistent with the observed optical galaxy luminosity func- 
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tion and thus does not exclude the existence of such objects. 
Some of the more luminous galaxies may have a clumpy dis- 
tribution of dust, which might allow for a larger fraction of 
their ionizing photons to escape as compared to other galax- 
ies. There could be an additional population of infrared se- 
lected galaxies that would not be included within the lumi- 
nosity function simulated here. Although these galaxies also 
might make some contribution to the ionizing background, 
a fairly small fraction of ultraviolet photons are thought to 
escape from them (Goldader et al. 2002). 

Assuming the conversion between the intensity at the 
Lyman limit and the one-sided flux seen by galaxies, as de- 
fined in Tumlinson et al. (1999), the first simulation, as 
illustrated in the figures, gives rise to a number of Ly- 
man limit absorbers per unit redshift {dN/dz)o,LL ~ 3.67 
when the number density of galaxies is adjusted to pro- 
duce {dN/dz)Q = 24.3 for forest absorbers as in Bahcall 
et al. (1996). The number density of galaxies used here is 
found to remain consistent with observed galaxy luminosity 
functions as discussed in Linder (1998). The observed values 
for (dN / dz)o^LL tend to be ~ 1 or lower (Lanzetta, Wolfe, 
& Turnshek 1995b; Storrie Lombardi et al. 1994; Stengler- 
Larrea et al. 1995). Adjusting the switching radius does not 
substantially change the number of Lyman limit absorbers, 
as these absorbers still arise largely in the exponential parts 
of the discs. However, if the scale length ratio hni/hB is de- 
creased to 1.2 then a more realistic number of Lyman limit 
absorbers arises as seen in Table |5| In this case an extra 
population of weaker absorbers would be needed in order 
to produce the observed {dN/dz)o. However most Lya ab- 
sorbers, including even the weakest ones, are found to trace 
the large scale galaxy distribution, so that they are likely 
to arise in gas which is about as close to galaxies as that 
in the first simulation. However this absorbing gas could be 
higher above the planes of galaxy discs and thus be even 
more highly ionized, as in the fifth simulation shown in Ta- 
ble 121 Alternatively the additional absorbers could behave 
more like randomly distributed points, in which case they 
would have a distribution of ionization rates similar to that 
for the first simulation, as illustrated in Fig. Q 



Another possibility is that Lyman limit absorbers arise 
only around luminous, high surface brightness galaxies. Ob- 
servers have also questioned whether galaxies which are low 
in luminosity and/or surface brightness give rise to Lyman 
limit absorption (McLin, Giroux, & Stocke 1998; Steidel 
1995; Steidel et al. 1997; Bergeron & Boisse 1991), although 
such objects are being found to give rise to stronger damped 
Lya absorption (Cohen 2001; Turnshek et al. 2000; Bowen, 
Tripp, & Jenkins 2001). Some possible explanations include 
gas being blown out more easily from less massive galaxies as 
suggested by McLin et al. (1998), differences in the column 
density proflles in the outer parts of LSB galaxies compared 
to what is simulated here, a slightly steeper neutral column 
density distribution compared to the N^]'^ assumed here, 
and/or a galaxy surface brightness distribution which allows 
for fewer moderate to large sized LSB galaxies compared to 
the flat surface brightness distribution simulated here. Other 
possible scenarios and implications for the nature of Lya ab- 
sorbers will be discussed further in a future paper. 
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Figure 1. A distribution of values for the direction and frequency 
averaged ionization rate is shown for simulated Lyo absorbers 
with N[jj > 10^^ cm~^ (solid line) and for randomly chosen 
points within the simulation (dashed line). Also shown (as *) 
is the assumed minimum background value due to quasars based 
upon Dave et al. (1999) and Haardt & Madau (1996). The triangle 
indicates the value for z < 1 found by Scott et al. (2002) using 
proximity cfTcct measurements. 



3 IONIZATION RATE FLUCTUATIONS 

A further understanding of the fluctuations in the ionizing 
background intensity due to galaxies is important for nu- 
merous reasons. First of all, the fluctuations need to be un- 
derstood as a source of uncertainty in measuring the overall 
background intensity. Histograms of values for the ionization 
rate are shown in Fig. Q for the first simulation, and are 
shown both for Lya absorbers {Nhi > 10^^ cm~'^) and for 
random points in the box. The histograms are shown along 
with the z — quasar value from Dave et al. (1999), based 
upon Haardt & Madau (1996) spectra, which is the assumed 
minimal value for Cj". 

Absorbers see higher values somewhat more often 
than randomly chosen points, as absorbers arise on average 
closer to galaxies. Overall the ionization rate seen by Lya 
absorbers ai z — varies over about a factor of about two 
in this simulation. None of the absorbers or random points 
are exposed to ionization rates which are within about ten 
percent of the minimal assumed background, although ex- 
tinction from gas between the absorbers and their ionization 
sources is not modelled here, so that in reality there may be 
a few. Most often both the random points and absorbers 
have values which are about a factor of 1.4 larger than (^i,. 
Varying the spectral index of the contribution from these 
nearby galaxies to 2 would simply move this peak to 

^ ~ l.S^b- Points with the largest values of (" tend to arise 
at quite low impact parameters from the centres of galaxies, 
as will be seen in Fig. |1| 

Also shown in Fig. is the measurement from Scott et 
al. (2002) for z < 1. This measurement has error bars which 
are large compared to the range of values shown in this 
plot, and the lower error bar would be at a slightly higher 
value than the peak of the plotted histogram, although an 
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Table 1. Ionization Rates for Varied Limiting Absorber Column 
Densities 



-13.5 -13.4 -13.3 -13.2 -13.1 

log(C/s"') 

Figure 2. A distribution of values for the ionization rate <^ is 
shown for simulated Ly« absorbers with varying limiting neutral 
column densities. Values of ( are binned for absorbers having 
;m ~2 (solid line), > 10^^ cm (long-dashed line), 
(dashed line), > lO^^ cm (dot-dashed line), and 
' (dotted line). Stronger absorbers arise on average 
closer to luminous galaxies, so that their f values are typically 
higher. 
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extra contribution from star-forming galaxies would move 
the calculated peak within the measurement error bars. 
While the galaxies are simulated here at z ~ 0, the ioniza- 
tion rate is seen to evolve only by < 0.2 orders of magnitude 
between the Scott et al. (2002) measurements at z < 1 and 
z > 1. According to the redshift evolution model they fit, 
where the background intensity evolves as (1 + z)"'^^'' , there 
would be even less difference expected between their mea- 
surement for z < 1 and what would be expected at 2; ~ 0. 
This could mean that more ionizing radiation escapes from 
galaxies than is assumed here, although it is more likely 
that this power law model does not describe the evolution 
in the ionizing background very well down to z ~ 0. How- 
ever, an additional uniform contribution could also come 
from galaxies at redshifts ^ 2, as the universe is optically 
thin to ionizing photons at these redshifts. Such a contri- 
bution would likely be ~ ("t, as calculated for star-forming 
galaxies (Giallongo et al. 1997; ShuU et al. 1999; Bianchi 
et al. 2001). Another possible contribution could come from 
gas in the intragroup medium (Maloney & Bland-Hawthorn 
1999; 2001). On the other hand, the other previously men- 
tioned low redshift measurements also tend to be consistent 
with ^ values closer to the assumed ("j,, so the value at 2: = 
need not be as large as the Scott et al. (2002) measurement. 

While the background intensity detected using the prox- 
imity effect is likely to be the most common value in a peaked 
distribution like the one found here, other values could be 
seen, for example, when making measurements of the ion- 
ization rate around a galax;y in an unusual environment. 

Distributions of ( values for absorbers are shown in 
Fig. I^l where the limiting neutral column density for ab- 
sorbers is varied. It can be seen that more high ^ values 
arise for stronger absorbers, as these arise on average closer 
to galaxies. The average logC value increases only slightly 
with limiting Nhi however, but the distributions become 
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For Lya absorbers with varying limiting neutral column densi- 
ties, Nj{j „^i„, the mean values for logf and standard deviations 
are shown for the first simulation as also illustrated in Fig.|2] Ab- 
sorbers with higher Nuj appear to have more higher f values in 
the Figure. Here it can be seen that the mean log ^ values change 
only slightly with limiting Njji, but the a values become much 
larger, indicating more substantial tails of high C, values. 




Figure 3. Values of the neutral column density are plotted ver- 
sus ionization rate for simulated Lya absorbers. Absorbers are 
shown as arising from 15,000 lines of sight in the region where 
Nuj < lOi^ cm~2 and ( > 10~13'3 s~l, and ten percent of these 
simulated points are shown in the region where Nhi < lOi^ cm~2 
and f < 10~i3-3 s~l which would otherwise be further saturated. 
Absorbers are shown as arising from 60,000 lines of sight through 
the box where N^i > lOi^ cm~2 Fairly high values of can be 
seen for any value of Nhi- Few points arise with Nhi ~ 10^^ 
cm~2 because absorbers change from neutral to highly ionized 
at a slightly larger Nhi value, as seen in galaxy disc ionization 
edges. 



much less strongly peaked due to larger tails of high val- 
ues, as shown in Tabled 

Values of ( are plotted versus neutral column density 
(Nhi) for absorbers in Fig.|3 Different numbers of simulated 
points are plotted in different regions of the plot in order to 
reduce saturation for low values and to show more detail 
for absorbers with high Nhi- Again the minimum assumed 
^ = produces a cutoff on the left side of the plot. Few 
points are seen with Nhi ~ lOi* to 10^^ , which is related to 
the galaxy disc ionization edges where Nhi falls ofi^ quickly 
with increasing radius around these values. It can be seen 
that absorbers with a wide range of Nhi can be exposed 
to a wide range of ionization rates. Even more absorbers 
with a wide range of higher values would be seen for a 
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Table 2. A Summary of the Simulations 



Description 


(dN/dz)o 


(dN/dz)o^LL 


log{C/s-i) 




1: As in text 


25.93 


3.91 


-13.38 


0.06 


2: T = 2.8 


27.43 


4.01 


-13.43 


0.05 


3: T = 3.5 


28.85 


4.09 


-13.47 


0.03 


4: T = 1.5 


21.72 


3.69 


-13.23 


0.10 


5: warp= 25° 


24.47 


3.69 


-13.36 


0.07 


6: h2i/hB = 1.2 


14.06 


0.90 


-13.44 


0.04 



The first simulation is as described in the text and illustrated 
in the Figures. Shown for each simulation in the table are the 
parameters which were varied (where each simulation is otherwise 
the same as the first simulation) , and the numbers of absorbers > 
-|^q14.3 (,ijj-2 g^j^j Lyman limit absorbers arising per unit redshift 
in a simulation with 12590 galaxies in a 28.9 Mpc cube, and the 
mean values for the logarithm of the ionization rate for Lya 
absorbers, and the standard deviation for the distribution of log C,. 

wide range of Nhi if more ionizing photons were allowed to 
escape from galaxies in the simulation. This happens even 
though one might expect absorbers exposed to large C, to be 
ionized away or seen only at lower values of Nhi- High Nhi 
values with high C, still arise, as the galaxies are simulated 
with a full range of disc inclinations. 

Occasionally voids are reported in the Lya forest 
(Crotts 1987; Dobrzycki & Bechtold 1991; Cristiani et 
al. 1995), and there is the possibility that variations in the 
ionizing background might contribute to such voids. How- 
ever at this time numbers of voids beyond what would be 
expected for a random absorber population have not yet 
been detected for voids smaller than the box size simulated 
here (Williger 2002). 



4 MODEL PARAMETER VARIATIONS 

Several model parameters were varied in further simulations 
as described in Table |21 The number of Lyman alpha ab- 
sorbers having Nhi > lO^'*'^ cm~^, the number of Lyman 
limit systems, and the mean and standard deviation a for 
log ( are shown for each simulation. In each case the mode 
for the distribution of log ( is very close in value to the mean. 
In cases where a is relatively large, there tends to be a more 
substantial tail of points having high ( values. The optical 
depth r was varied in order to explore the uncertainty range 
in the fraction of ionizing photons escaping from galaxies. 
The value of r = 2.8, preferred for our Galaxy by Bland- 
Hawthorn (1998), corresponds to a direction-averaged es- 
cape fraction of ~ 2 percent of the ionizing photons (Bland- 
Hawthorn & Maloney 2001). In the first simulation, it is 
assumed that r is dependent upon galaxy central surface 
brightness as discussed above, while in the second simula- 
tion we assume r = 2.8 for all galaxies. In the third simula- 
tion r — 3.5 is used for all galaxies in order to produce an 
escape fraction of ~ 1 percent, while r = 1.5 is used for all 
galaxies in the fourth simulation, giving an escape fraction 
~ 10 percent. A similar distribution for arises when the 
value T = 2.8 is assumed for all galaxies as compared to the 
first simulation. The mean and mode values of log ( corre- 
spond to 1.38(^6 in the first simulation, 1.23(^6 in the second 
simulation, Lll("5 where r — 3.5, and 1.95("(, where r — 1.5. 



Thus the uncertainty in the fraction of ionizing photons es- 
caping from galaxies means that normal galaxies could con- 
tribute between ten percent of the quasar background and 
an amount about equal to the quasar background. 

The typical disc warping angle, which is also rather un- 
certain, was also varied as shown in the fifth simulation in 
Table|21 Bland-Hawthorn (1998) has suggested, for example, 
that there may be a selection bias against detecting highly 
warped discs in HI because they become more highly ion- 
ized by the stars within the galaxy. When the disc warping 
angle is increased to 25° then a slightly higher mean log 
is seen, as more ionizing photons escape higher above the 
plane of a disc. However the absorbers even in this case are 
not generally far above the plane of the disc, as the disc is 
only warped beyond two HI scale lengths. 

The distributions in ionization rates seen in the simu- 
lations here are generally quite strongly peaked, although 
there may be more absorbers with high if there is sub- 
stantial absorbing gas above the planes of galaxy discs. The 
contribution to the ionization rate varies by a factor of 45 
for the Bland-Hawthorn (1998) model between polar angles 
of 0° and 80° at some distance from a galaxy with r = 2.8. 
An additional simulation was done where absorbers arise in 
galaxy haloes with column density profiles obeying equa- 
tion (24) from Chen et al. (1998). While the absorbing gas 
is not modeled in this case, the distribution of ionization 
parameters is found to be similar to that seen in the fifth 
simulation (where the warping angle = 25°). Even in this 
case, however, or for randomly chosen points, or in any case 
where absorbers arise in random directions relative to galax- 
ies, few of the points or absorbers will arise close to the 
galactic poles. 

The sixth simulation was done with a reduced /i2i/^s = 
1.2 in order to illustrate a scenario with a more realistic 
number of Lyman limit absorbers. In this case the mean log C, 
is slightly lower compared to the first simulation because 
the simulated absorbers arise typically somewhat closer to 
the centres of galaxies within the warped outer discs, so 
that they arise closer to the planes of the discs. In order 
to make such a scenario realistic, however, an additional 
population of low column density absorbers would be needed 
to account for the observed {dN/dz)o as in Bahcall et al. 
(1996). The additional absorbers could either arise far from 
luminous galaxies and thus behave more like the random 
points illustrated in Fig.0 or they could arise higher above 
galactic planes where they could be more highly ionized. 

A simulation was also run at redshift one by decreas- 
ing the box size by a factor of (1 + z) and assuming ^5 = 
3.33 X 10~^^ s~^, again based upon the calculations of Dave 
et al. (1999) and Haardt & Madau (1996) at z = 1. Any 
inconsistency might indicate some galaxy luminosity evolu- 
tion and/or require more photons to escape from galaxies at 
z — 1, although the ^ distribution was found to be peaked 
at 3.7 X 10"^^ s"\ which is between the Scott et al. (2002) 
measurements for z < 1 and z > 1. 



5 RELATIONSHIP TO GALAXIES 

Understanding how the ionizing background varies with 
galactic environment will be important for understanding 
what kinds of objects give rise to Lya absorption, and in 
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Figure 4. The average impact parameter to the nearest four 
galaxies having Mb < —16 is plotted versus ionization rate for 
simulated Lyo absorbers with Nhi > lO'^^ cm~-^. Circles indicate 
points where the nearest galaxy having Mb < —16 is within 
100 kpc, and triangles indicate such a galaxy within 200 kpc. 
Eight percent of the simulated points with the nearest galaxy 
> 200 kpc and fifteen percent of the other points are shown in the 
region where f < 10~^^'^ s~^ where the points would otherwise 
be further saturated. Absorbers with the largest ( values tend to 
be in environments which are rich in nearby galaxies. Absorbers 
with lower values of ( arise at a wide range of averaged impact 
parameters. 



Figure 5. Absolute magnitudes are plotted for galaxies where 
Lya absorption (> lO^'^ cm~^) is known to arise in the simula- 
tion versus ionization rate for the absorber. Ten percent of the 
simulated points are shown for logf < —13.3 where the figure 
would otherwise be further saturated. Horizontal lines are seen 
because numerous absorbers can arise close to a particular sim- 
ulated luminous galaxy. More high ^ absorbers are seen around 
luminous galaxies simply because these galaxies are assumed to 
have more absorbing gas around them so that more absorbers 
with any f value arise. Fairly high values of can arise around 
galaxies with any Mb value. 



which environments Lya absorption is most likely to arise. 
Tripp et al. (1998) suggest that no absorption is found near 
a galaxy cluster due to increased ionization. Fluctuations in 
the ionizing background may also allow for substantial vari- 
ations in the environments affecting galaxy formation pro- 
cesses. For example, it has been suggested that dwarf galax- 
ies may form less easily in an intense radiation field (TuUy et 
al. 2002; Efstathiou 1992; Quinn, Katz, & Efstathiou 1996; 
but see the discussion in Sabatini et al. 2003). 

In order to attempt to parametrize the galaxy clustering 
environment, in Fig.QlC values are plotted (again for the first 
simulation) versus the average distance to the nearest four 
galaxies having Mb < —16, as observers are generally un- 
able to detect less luminous galaxies around even the nearest 
absorbers. Again here and in the next two figures, different 
numbers of simulated points are plotted for log^ < —13.3 
values in order to reduce the saturation in the plot. The near- 
est four galaxies are used because the ionizing background 
is likely to be somewhat higher even in a group environment 
in addition to being higher in a rich cluster. The shapes of 
the plotted points give more information about the near- 
est single galaxy with Mb < —16. The ( values are shown 
for absorbers with Nhi > 10^^ cm~'^. Most of the points 
with C ~ Cb &re far from luminous galaxies, while those with 
higher arise more often closer to luminous galaxies. How- 
ever, selection effects against LSB galaxies may make this 
correlation less clearly visible to an observer. Some higher 
values can arise even far from luminous galaxies, however, 
as it was assumed that some ionizing radiation escapes even 
from dwarf galaxies which are not assumed to be strongly 
clustered. Still the presence of any luminous galaxy may 
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Figure 6. Central surface brightnesses are plotted for galaxies 
where Lya absorption (> lO^'^ cm^^) arises nearby versus ion- 
ization rate for the absorber. Ten percent of the simulated points 
are shown for logC^ < —13.3 where the figure would otherwise be 
further saturated. High f values can arise around galaxies with a 
wide range of central surface brightness values. 



have a more important effect on the ionization rate rather 
than the overall clustering environment, as large values 
tend to arise more often when a luminous galaxy is within 
200 kpc. 

In Figs. |S] and |H] ^ values are plotted versus galaxy lu- 
minosity and central surface brightness for each absorber, 
where their associated galaxies are known from the sim- 
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Figure 7. The average ionization rate, shown in units of the 
assumed quasar background f;,, is shown here for absorbers 
> 10^^ cm~-^ versus impact parameter to the nearest galaxy with 
Mg > —16. A solid line is shown where the nearest such galaxy 
has Mb < -19, a dashed line for -19 < Mb < -18, dot-dashed 
for -18 < Mb < -17, and dotted for Mb < -17. It can be seen 
that the average f value for absorbers decreases with increasing 
impact parameter from a luminous galaxy, and the decrease ap- 
pears to be steepest for the most luminous galaxies. Note that 
the error bars based upon the standard deviation in simulated 
values would be smallest for impact parameter values in the cen- 
tre of the plot where C/Cfc ~ 1-4, as there are fewer data points 
for high and low f values. 



ulation. (Note that an observer might identify a different 
galaxy as associated with an absorber, as discussed in Lin- 
der (2000).) Points are seen to lie on horizontal lines in ei- 
ther plot, as galaxies which are luminous or moderately low 
in surface brightness are assumed to have large absorption 
cross sections for their surrounding gas and give rise to nu- 
merous absorbers within the simulation. It can be seen in 
either plot that absorbers arising around galaxies with a 
wide range of properties are exposed to a wide range of ion- 
ization rates. Thus variations in ^ happen around particular 
galaxies, although these galaxies can have a wide range of 
properties. High ^ values will be seen most often in locations 
where absorbers arise most often, such as those close to lu- 
minous galaxies. Yet many more faint galaxies exist, where 
some absorbers with high values can also be seen, and no 
evidence is seen for a variation in the average values with 
galaxy luminosity or surface brightness. 

The ionization rate tends to be higher close to galaxies 
and in regions of higher galaxy density where Lya absorbers 
often arise, but how much is the intergalactic medium af- 
fected on average by ionizing radiation from galaxies? The 
lowest ^ values tend to be seen when looking as far as possi- 
ble from a luminous galaxy, as can be seen in Fig.|7| A fall-off 
can be seen in the average value with impact parameter 
from a galaxy with Mb < —16, and such a fall-off appears to 
be steeper for more luminous galaxies. Again this plot may 
be affected by luminous galaxies having more absorbers with 
high ^ values around them simply because more absorbing 
gas is assumed to be located around luminous galaxies. 

A less biased view of the ionization of the intergalactic 
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Figure 8. The average ionization rate, shown in units of the 
assumed quasar background (^t,, is shown here at randomly cho- 
sen points versus impact parameter to the nearest galaxy with 
Mb > —16. A solid line is shown where the nearest such galaxy 
has Mb < -19, a dashed line for -19 < Mb < -18, dot-dashed 
for -18 < Mb < -17, and dotted for Mb < -17. Compared 
to the previous figure, it can be seen even more clearly here that 
the average value for absorbers decreases with increasing im- 
pact parameter from a luminous galaxy, and the decrease is the 
steepest for the most luminous galaxies. 

medium can be seen from looking at a similar plot of the av- 
erage C value versus impact parameter to the nearest galaxy 
for random points rather than absorbers, as shown in Fig.|H] 
While it becomes even more difflcult here to simulate points 
which are very close to galaxies, it can be seen even more 
clearly that the intergalactic medium is more highly ionized 
at points which are closer to more luminous galaxies. While 
Fig.Qis a prediction of what observers might see if they are 
able to measure for numerous absorbers. Fig. |S]is a bet- 
ter representation of how the ionization of the intergalactic 
medium could be simulated. 



6 CONCLUSIONS 

Normal galaxies are likely to contribute at least thirty to 
forty percent of what quasars do to the ionizing background 
of Lyman continuum photons at zero redshift, assuming that 
ionizing photons escape from other galaxies in an analogous 
manner to the Bland-Hawthorn (1998) model of our Galaxy 
where ~ 2 percent of ionizing photons escape. Allowing for 
some uncertainty in this direction-averaged ionizing photon 
escape fraction, assuming that between ~ 1 and ~ 10 per- 
cent of ionizing photons escape means that the contribution 
to the ionizing background from normal galaxies could be 
between ~ 10 and ~ 100 percent of the assumed quasar 
contribution. This ultraviolet background is important for 
ionizing gas surrounding galaxies and within the intergalac- 
tic medium. This gas gives rise to Lya absorption and makes 
an uncertain contribution to the baryon content in the local 
universe due to uncertainties in ionization intensities and 
mechanisms. 

Distant quasars at somewhat higher redshifts will ion- 
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ize the low redshift universe in a relatively uniform man- 
ner, but ionizing radiation escaping from normal galaxies at 
low redshift will result in local fluctuations in the ionizing 
background. Fluctuations have been found in the ionization 
rate of gas around sinmlatcd galaxies, which will give rise to 
variations in the neutral gas fractions in Lya absorbers with 
a wide range of neutral column densities. A wide range of 
ionization rates are found to arise close to galaxies having 
a wide range of properties, but normal galaxies also play 
an important role in ionizing the more distant intergalac- 
tic medium. Ionization rates for absorbers are found to be 
about twice as high as the quasar background on average 
when looking at ^ 200 kpc from a luminous galaxy, about 
1.4 times the quasar background when looking at ~ 1 to 
2 Mpc from a luminous galaxy, and only as little as ~ 1.1 
times the quasar background when looking as far as possi- 
ble from luminous galaxies. Luminous star forming galaxies, 
which may contain rather clumpy dust and thus have more 
complex extinction behaviour than was modelled here, will 
also give rise to both local and larger scale fluctuations in 
the ionizing background. 

Fluctuations in the ionizing background may have im- 
plications for the formation and evolution, and our ability 
to detect, smaller objects located near luminous galaxies, 
such as dwarf galaxies and high velocity clouds. Such fluc- 
tuations will also give rise to variations in galaxy disc ion- 
ization edges, which will have implications for the column 
density distribution of Lyman limit absorbers. Fluctuations 
in the ionizing background may also have some implications 
for the nature of Lya absorbers and their relationship to 
galaxies. If galaxies play am important role in ionizing the 
gas around them, gas which must make some contribution 
to Lya absorption, then further constraints may be made 
on models for galEixies giving rise to absorption. Gas which 
is too close to a luminous galaxy, and in particular that lo- 
cated far above the plane of the disc, will be exposed to 
more ionizing radiation than absorbing gas which is ion- 
ized largely by a quasar background, possibly reducing the 
number density of absorbers that will arise in such environ- 
ments. Variations in the fractions of ionizing radiation that 
escape from galaxies with various properties may be impor- 
tant for determining what kinds of galaxies can give rise to 
absorption. A further understanding of the distribution of 
ionization rates will enable us to learn more about chemical 
abundances using metal absorption line systems. 

The fluctuations in the ionizing background arc not 
seen to be substantial in the simulations here however, al- 
though most absorbers are assumed to arise fairly close to 
the planes of galaxy discs. Larger variations in the photoion- 
ization rates would only be seen if substantial amounts of ab- 
sorbing gas are concentrated above galactic poles, although 
in this case the gas might include components which are 
ejected from the galaxies so that coUisional ionization pro- 
cesses might also be important. 
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APPENDIX A: RELATIONSHIP BETWEEN 
THE COLUMN DENSITY PROFILE IN A DISC 
AND THE COLUMN DENSITY DISTRIBUTION 

Suppose absorbers arise in the outer parts of galaxy discs, 
where the column density Nhi in each disc falls off with 
radius r as a power law where Nhi oc r~^. The neu- 
tral column density distribution dAf/dNni resulting from 
these absorbers in then dV /dNni = dN jdr x dr/dNui = 
2nr X dr/dNni- Since dNni/dr oc r~^^~^^^ from the as- 



sumed column density profile then dM /dNni 



r.{P+2) 



[Njjj'^Y^^'^K Thus when the neutral column density dis- 



-l/p\(p+2) 

II ) 

tribution in each outer disc falls off as a power law with 
exponent —p, then the resulting neutral column density dis- 
tribution is also a power law having exponent — e, where 
dU/dNui oc N-] oc N~^^+'^^''' and p = 2/(e - 1). 



APPENDIX B: POLAR ANGLE CALCULATION 

At a point (x, y, z) in space where the ionization rate is 
calculated, each galaxy i contributes emission which is seen 
at angle 0i 7r/2) from the galaxy's pole. Each galaxy has 
a randomly simulated inclination 6 between the disc plane 
and the a-axis, thus chosen from a uniform distribution in 
cos 6', and a random orientation (j>, where ^ (f> <2ii. 

The angle between the galaxy's rotation axis given by 
the vector v and the line given by the vector d = (x — Xi,y — 
yi,z — Zi) equals 
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This angle is always in the interval [0, 7r/2]. Here ||d|| = 

yfdl+dl + dl. The vector V is produced by rotating 
the vector (0,0,1) first by n/2 — 9 in the xz-plane, 
which gives the vector (cosS, 0, sinS), and then rotat- 
ing by angle in the a;j/-plane, which gives the vector 
(cos ^ cos 0, — cos 6 sin 0, sin 0) . Hence ||v|| = 1. 
The polar angle is thus 



\dx cos 9 cos (f) — dy cos 9sm(l> + dz sin 9\ 

(diTdfTdlF^ 



Gi = cos ^ 
where 

dy = y-yi 
dz = z — Zi 

for galaxy i centred at {xi,yi,Zi). 
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